EFFECTIVE MACROSCOPIC DYNAMICS OF STOCHASTIC 
PARTIAL DIFFERENTIAL EQUATIONS IN PERFORATED 

DOMAINS 

WEI WANG, DAOMIN CAO & JINQIAO DUAN 

Abstract. An effective macroscopic model for a stochastic microscopic 
system is derived. The original microscopic system is modeled by a sto- 
chastic partial differential equation defined on a domain perforated with 
small holes or heterogeneities. The homogenized effective model is still a 
stochastic partial differential equation but defined on a unified domain with- 
out holes. The solutions of the microscopic model is shown to converge to 
those of the effective macroscopic model in probability distribution, as the 
size of holes diminishes to zero. Moreover, the long time effectivity of the 
macroscopic system in the sense of convergence in probability distribution, 
and the effectivity of the macroscopic system in the sense of convergence in 
energy are also proved. 



1. Introduction 

In recent years there has been explosive growth of activities in multiscale 
modeling of complex phenomena in many areas including material science, cli- 
mate dynamics, chemistry and biology [151 [31]. Stochastic partial differential 
equations (SPDEs or stochastic PDEs) — evolutionary equations containing 
noises — arise naturally as mathematical models of multiscale systems under 
random influences. In fact the need to include stochastic effects in mathemat- 
ical modelling of realistic physical behavior has become widely recognized in, 
for example, condensed matter physics, climate and geophysical sciences, and 
materials sciences. But implementing this idea poses some challenges both in 
theory and for computation [T71 1 
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This paper is devoted to the effective macroscopic dynamics of microscopic 
systems modeled by parabolic SPDEs in perforated media which exhibit small- 
scale heterogeneities. One example of such microscopic systems of interest 
is composite materials with microscopic heterogeneities under the impact of 
external random fluctuations. 

The heterogeneity scale is taken to be much smaller than the macroscopic 
scale, which is equivalent, here, to assuming that the heterogeneities are evenly 
distributed. From a mathematical point of view, one can assume that micro- 
scopic heterogeneities (holes) are periodically placed in the media. This peri- 
odicity can be represented by a small positive parameter e (i.e., the period). 
In fact we work on the space-time cylinder D e x (0, T), with T > 0, and D t 
being the spatial domain obtained by removing a number N € of holes, of size 
e, periodically distributed, from a fixed domain D. When taking e — > 0, the 
holes covering D are smaller and smaller and their numbers goes to oo. This 
signifies that the heterogeneities are finer and finer. 

There are lots of work on the homogenization problem for the deterministic 
systems defined in such perforated domain or other heterogeneous media, see 
for example [01 [2H E51 [2H ED] for heat transfer in a composite material, [21 El EI] 
for the wave propagation in a composite material and [211 123] for the fluid flow 
in a porous media. For an introduction see [HI [27] . 

Recently there are also works on homogenization in the random context; see 
[T9l [22| [26] for general random coefficients, and [51 [351 EE] for randomly perfo- 
rated domains. And also see a survey book about the homogenization results 
in a random context [T5]. A basic assumption in these texts is the ergodic 
hypotheses on the coefficients for the passing of the limit of e — > 0. Note that 
the microscopic models in these works are partial differential equations with 
random coefficients, so-called random partial differential equations (random 
PDEs), instead of stochastic PDEs. 

In the present paper, the microscopic model is a SPDE defined in a per- 
forated domain. Homogenization techniques are employed to derive an effec- 
tive, simplified, macroscopic model. Homogenization is a formal mathematic 
procedure for deriving macroscopic models from microscopic systems. It has 
been applied to a variety of problems including composite materials modelling, 
porous media and climate modelling; see [HI EDI ESI EZ] • Homogenization pro- 
vides effective macroscopic behavior of the system with microscopical hetero- 
geneities for which direct numerical simulations are usually too expensive. 

We consider a spatially extended system where stochastic effects are taken 
into account in the model equation, defined on a deterministic domain but 
perforated with small scale holes. Specifically, we study a class of stochastic 
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partial differential equations driven by white noise on a perforated domain in 
the following form 

du e {t) = (A e u e + F e (x, i))dt + G e (x, i)dW{t), < t < T, e > 0, 

which will be described in detail in the next section. For the general theory of 
SPDEs we refer to p2]. The goal here is to derive the homogenized equation 
(effective equation), which is still a stochastic partial differential equation, for 
([T]) by the homogenization techniques in the sense of probability. 

Homogenization theory has been developed for deterministic systems, and 
compactness discussion for the solutions {u e } e in some function space is a key 
step in the homogenized approach [!9J. However, due to the appearance of the 
stochastic term in the above microscopic system considered in this paper, such 
compact result does not hold for this stochastic system. Fortunately the com- 
pactness in the sense of probability, that is, the tightness of the distributions for 
{u e }, still holds. So one appropriate way is to homogenize the stochastic sys- 
tem in the sense of probability. The goal in this paper is to derive an effective 
macroscopic equation for the above microscopic system, by homogenization in 
the sense of probability. It is shown that the solution u e of the microscopic 
or heterogeneous system converges to that of the macroscopic or homogenized 
system as e j in probability. That also implies the distributions of {u e } t 
weakly converge, in some appropriate space, to the distribution of a stochastic 
process which solves the macroscopic effective equation. Moreover, the long 
time effectivity of the homogenized macroscopic system is demonstrated, that 
is, the solution u e (t) is shown to converge to the stationary solution of the 
homogenized equation as t — > oo and e j, in the sense of probability distri- 
bution. Furthermore, the effectivity of the macroscopic system in the sense of 
convergence in energy is also shown. 

In our approach, one difficulty is that the spatial domain is changing as 
e — ► 0. To overcome this we use the extension operator introduced in [8] and 
introduce a new probability space depending on a parameter in which the so- 
lution is uniformly bounded. One novelty here is that the original microscopic 
model is a stochastic PDE, instead of a random PDE as studied by others, 

e.g., HSl|2BllZ|. 

This paper is organized as follows. The problem formulation is stated in 
§|2j Section [3] is devoted to basic properties of the microscopic system. The 
effective macroscopic equation is derived in §HJ The long time effectivity of the 
homogenized macroscopic system is considered in £j6l Finally, the effectivity of 
the macroscopic system in the sense of convergence in energy is shown in £j5l 
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Moreover, in the Appendix we present the explicit expression of the homoge- 
nization matrix. 



2. Problem formulation 

Let D be an open bounded set in R™, n > 2, with smooth boundary 3D 
and e > is a small parameter. Let Y = [0, Zi) x [0, Z 2 ) x • ■ ■ x [0, l n ) be 
a representative (cubic) cell in W 1 and S an open subset of Y with smooth 
boundary dS, such that S C Y. Write I = (h,h,-- - ,l n )- Define eS = 
{ey : y G S}. Denote by S e> k the translated image of eS by kl, k E Z n , 
kl = (k^i, k 2 l2, • • • , k n l n ). And let S e be the set of all the holes contained in 
D and D e = D\S e . Then D t is a periodically perforated domain with holes of 
the same size as period e. We assume that the holes do not intersect with the 
boundary dD, which implies that dD e = 3D U dS e . See Fig. 1 for the case 
n = 2. This assumption is for avoiding technicalities and the results of our 
paper will remain valid without this assumption; see pp. 

In the sequel we use the notations 



Y* = Y\S, d 



\Y* 



with \ Y\ and \Y*\ the Lebesgue measure of Y and Y* respectively. And denote 
by v the zero extension to the whole D for any function defined on D e : 



v on D e , 
on S f . 




Fig. 1: Geometric setup in M. 2 
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Now for T > fixed final time, we consider the following ltd type nonau- 
tonomous stochastic partial differential equation defined on the perforated do- 
main D e in M n . 

du e (x,t) = (div(A e (x)Vu e (x,t)) + f e (x,t)^dt + g e (t)dW(t) (2.1) 
in D e x (0,T), 

u e = on dD x (0,T), (2.2) 
^P- = on dS e x (0,T), (2.3) 
u e (0) = u° m D e , (2.4) 
where the matrix A £ is 



and 



9 • v-^ ( x \ 9 



az/4 V e / oxj 

ij j 

with n the exterior unit normal vector on the boundary dD e . 
We make the following assumptions on the coefficients: 

(1) Oy e L°°(M n ), z,j = l,...,n; 

(2) S"j- =1 a>ij£i£j > « X^r=i f° r £ e ^™ an d a a positive constant; 

(3) ay are V-periodic. 

Furthermore we assume that 

f e eL 2 (D t x[0,T]) (2.5) 
and for < t < T, g e (t) is a linear operator from i 2 to L 2 (D e ) defined as 

oo 

g e (t)k = gl(x, t)ki, k = (k u k 2 , ■ ■ ■ ) G £ 2 

i=l 

where g\(x, t) e L 2 (D t x [0, T]), i = 1, 2, • • • , are measurable functions with 

oo 

El^0M)lW)<^, te[o,T) (2.6) 

i=l 

for some positive constant Ct independent of e. In ( 12. ip . W(t) = (Wi(t), W2(t), ■ 
is a Wiener process in C 2 with covariance operator Q = Idp. and {Wi(t) : z = 
1, 2, • • • } are mutually independent real valued standard Wiener processes on a 



6 W. WANG, D. CAO & J. DUAN 

complete probability space (fi, J-,¥) with a canonical filtration (J-t)t>a- Then 

oo 

\9M# = l£\9fr,t)\b(D.)<CT, te[0,T]. (2.7) 
i=i 

Here is the space of Hilbert-Schmit operators [TJl [TO] . Denote by E the 
expectation operator with respect to P. 

The following compactness result [20] will be used in our approach. Let 
X C y C Z be three reflective Banach spaces and X C y with compact and 
dense embedding. Define Banach space 

G = {v : v e L 2 (0,T;X),^ E L\0,T;Z)} 

at 

with norm 




Lemma 2.1. If B is bounded in G, then it is precompact in L 2 (0,T;y). 

Let S be a Banach space and S' be the strong dual space of S. We recall the 
definitions and some properties of weak convergence and weak* convergence 

m- 

Definition 2.2. A sequence {s n } in S is said to converge weakly to s G S if 
W G S' , 

lim (s', s n ) S ',s = 0', s) S ',s 

n— too 

which is written as s n — s weakly in S. Note that (s', s) denotes the value of 
the continuous linear functional s' at the point s. 

Lemma 2.3. (Eberlein-Shmulyan) Assume that S is reflexive and let {s n } 
be a bounded sequence in S. Then there exists a subsequence {s Hk } and s G 
S such that s nk — ^ s weakly in S as k — > oo. If all the weak convergent 
subsequence of {s n } has the same limit s, then the whole sequence {s n } weakly 
converges to s. 

Definition 2.4. A sequence {s' n } in S' is said to converge weakly* to s' G S' 
ifVs G S, 

lim (s' n , s) S ',s = (V, s) S ',s 
which is written as s' n — ^ s' weakly* in S' . 

Lemma 2.5. Assume that the dual space S' is reflexive and let {s' n } be a 
bounded sequence in S' . Then there exists a subsequence {s' n } and s' G S 1 
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such that s' nk — ^ s' weakly* in S' as k — > oo. If all the weakly* conver- 
gent subsequence of {s' n } has the same limit s' , then the whole sequence {s' n } 
waekly* converges to s' . 

We also use the following definition of the weak convergence of the Borel 
probability measures on S, for more we refer to [H]. 

Definition 2.6. Let {fi e }e be a family of Borel probability measures on the 
Banach space S. We say \x e weakly converges to a Borel measure fi on S if 

/ hd[i e — > / hdfi, as t [ 0, 
Js Js 

for any h G C&(«S) ; the space of bounded continuous functions on S. 

In the following, for a fixed T > 0, we always denote by Ct a constant 
independent of e. 



3. Basic properties of the microscopic model 



In this section we will present some estimates of the solutions of microscopic 
model (12.11) . useful for the tightness result of the distributions of solution 
processes in some appropriate space. 

Let H = L 2 (D) and H e = L 2 (D e ). Define the following space 

V e = {ueH\D t ),u\ dD = 0} 

provided with the norm 

'x\ dv 



\v\v. = \^aM 



/ dxj 



i=l 



T,H e 



This norm is equivalent to the usual H l (D e )-norm, with an embedding con- 
stant independent of e, due to the assumptions on in the last section. Here 
Q) n denotes the direct sum of the Hilbert spaces with usual direct sum norm. 
Let 

( dv 
V(A e ) = <veV e : div(A e Vv) E H e and 







dS, 



and define operator A e v = div(A t Vv) for v G V(A e ). Then system (12.ip - (12.4p 
can be written as the following abstract stochastic evolutionary equation 

du e = (A e u e + f e )dt + g e dW, u e (0) = u°. (3.1) 

By the assumptions on a^, operator A e generates a strongly continuous 
semigroup S e (t) on H e . Solution of (13.11) can then be written in the mild sense 

u e (t) = S e (t)u° e + [ S e (t-s)f e (s)ds+ [ S e (t-s)g e (s)dW(s) (3.2) 
Jo Jo 
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And the variational formulation is 



(du e (t), v) H -i v = ( — / A € (x)Vu e {x, t)Vv (x)dx + / f t (x,t)v(x)dx^dt + 



g e {x,t)v{x)dW{t), in V'{0,T), v G V e , (3.3) 



D, 



with ■u e (0,x) = u®(x). 



For the well-posedness of system (13.11) we have the following result. 

Theorem 3.1. (Global well-posedness of microscopic model) Assume 
l\2.5\) and {2. 7| ) hold. Let be a (jF , £>(-£f e )) -measurable random variable. 
Then system k3. 1\) has a unique mild solution u G L 2 (Q, C(0, T; H e )r\L 2 (0, T; V e )) , 
which is also a weak solution in the following sense 

(u e (t),v) St 

= (u° e ,v) He + / (A e u e (s),v) H Js + / (f e ,v) He ds+ / (g e dW, v) He (3.4) 
Jo Jo Jo 

for t G [0, 7) and v G V e . Moreover, if u° e is independent of W{t) with 
E\u Q t \ 2 Ht < oo, then 

E\u e {t)\k + E / \u e (s)\ 2 v ds < E\u° t \ 2 He +C T , for t G [0,T], (3.5) 
Jo 

E / \u e (s)\ 2 H ^ds < C T (E\u% e + 1), for t G [0,7]. (3.6) 
Jo 

If further assume that 

oo 

|V^ e (t)|Jo = £ |V Ae ^(t)| 2 e „ He < C r , /or f G [0,T] (3.7) 

i=l 

wii/i 6 K e and E\u®\y e < oo, then 

E\u e (t)\ 2 Ve + E [ \AMs)\ 2 H ds<E\u%+C T , forte[0,7}. (3.8) 
Jo 

Moreover, system ( Iff. II) is well-posed on [0, oo) if 

/ e GL 2 (0,oo;ff e ), &e£ 2 (0,oo;£?). (3.9) 
Proof. By the assumption (12 .7ft . we have 

oo 

u(t)\\ 2 c , = J2\9i(t^)\k<^- 

i=l 
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Then the classical result of [12] yields the local existence of u e . And applying 
the stochastic Fubini theorem, it is easy to verify the local mild solution is also 
a weak solution. 

Now we give the following a priori estimates which yields the existence of 
weak solution on [0, T] provide fl2.5jl and (12.71) hold. 
Applying ltd formula to \u e \ 2 , we derive 

d\u £ (t)\ 2 He -2(A e u e} u e ) He dt= 2(f e} u e ) He dt + 2(9edW,u e ) He + \g e \l Q dt (3.10) 
By the assumption on a^, we see that 

-(A e u e ,u e ) He > M u e\ 2 H e 

for some constant A > independent of e. Then integrating (13.101) with respect 
to t yields 

W{t)? He + / W\yds 
Jo 

L 2 (0,T;ff E ) + J (gedW,U e ) He ds + J Ig^gds. 

Taking expectation on both sides of the above inequality, we derive (13. 5p . 

In a similar way, application of ltd formula to \u t \ Vc = |V A E Ue\% n H e results 
in the relation 

dW{t)\\ + 2(A e u e ,A e u e ) Hc dt 
= -2(f e ,A e u e ) He dt - 2(g e dW,A e u e ) He + \V Ac g e \l Q dt. (3.11) 

Integrating both sides of (13. lip and by the Cauchy-Schwarz inequality, it is 
easily to have 

\ u e(t)\ v , + / \A e u e \ 2 He ds 
Jo 

< \Ue(0)\v t + \Mh(0,T;Hc) ~ 2 ^ (dedW, A e U e ) Hf d S + J \V A^e^ds, 

Then taking the expectation, we derive (13. 8p . By (13.31) and the property of the 
stochastic integral we easily have (13. 6p . 

Thus, by the above estimates, the solution can be extended to [0, oo) if (13. 9p 
hold. The proof is complete. □ 

We recall a probability concept. Let z be a random variable taking values 
in a Banach space S, namely, z : Q — > z. Denote by C(z) the distribution (or 
law) of z. In fact, C(z) is a Borel probability measure on S defined as [12] 

C(z)(A) = ¥{uj : z(u) e A}, 
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for every event (i.e., a Borel set) A in the Borel a— algebra B(S), which is the 
smallest a— algebra containing all open balls in S. 

As stated in £JU, for the SPDE (12. ip we aim at deriving an effective equa- 
tion in the sense of probability. A solution u e may be regarded as a random 
variable taking values in L 2 (0,T; H t ). So for a solution u e of (12. II) - (12. 4p de- 
fined on [0,T], we focus on the behavior of distribution of u e in L 2 (0, T; H € ) 
as e — > 0. For this purpose, the tightness [H] of distributions is needed. Note 
that the function space changes with e, which is a difficulty for obtaining the 
tightness of distributions. Thus we will treat {C(u e )} t> o as a family of distri- 
butions on L 2 (0,T; H) by extending u € to the whole domain D. Recall that 
the distribution (or law ) of u e is defined as: 

C(u e )(A) = F{u : u t (-,-,uj) G A} 

for Borel set A in L 2 (0, T; H e ). First we define an extension operator P e in the 
following lemmas. 

In the following we denote by h(X, y) the space of bounded linear operator 
from Banach space X to Banach space y. 

Lemma 3.2. There exists a bounded linear operator 

Q G h(H k (Y*) , H k (Y)) , k = 0,l, 

such that 

\VQv\® nL 2 {Y) < C\Vv\® nL 2 {Y *), v G H\Y*) 
for some constant C > 0. 

For the proof of Lemma 13.21 see [8] . 

We define an extension operator P e in terms of the above bounded linear 
operator Q in the following lemma. 

Lemma 3.3. There exists an extension operator 

P e G L(L 2 (0, T; H k (D e )), L 2 (0, T; H k (D))) , k = 0,1, 

such that for any v G H k (D t ) 

(1) P e v = v onD £ x (0,T) 

(2) \P e V \l 2 (0,T;H) < Ct\v \L 2 (0,T;H e ) 

(3) |VA £ (-Pe^)|L2(0,T;enL 2 (D)) < C*T | V^ E ^ | L 2 (0,T;ffi„L 2 (D e )) 

where Ct is a constant independent of e. 
Proof. For ip G H k (D t ), then 

<Pe{y) = ~<p{e y) 
e 
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belongs to H k {Y l *) with Y* the translation of Y* for some I £ R n . Define 

Q e ip{x) = e{Qip e ){-). (3.12) 
Now for up £ L 2 (0,T; H k (D e )), we define 

(p eV 0(*,t) = [4</>(v)](f) =e[g^,-)](f)- 

It is known [8] that the operator P e £ L(L 2 (0, T; H k (D e )), L 2 (0, T; H k (D))) , jfc 
= 0, 1, and satisfies the conditions (l)-(3) listed in the lemma. This completes 
the proof. □ 

Remark 3.4. In Lemma 2.1 of [8j ; the operator P t defined in L(L°°(0,T; 
H k (D t )), L°°(0, T; H k (D))) , k = 0, 1, coincides with the operator defined in 
Lemma lSTSi above. 

Remark 3.5. The estimates in Theorem \3. II for also hold for P e u e . In fact 
estimates ( Iff. 5\) and / Iff. 8\) are easily derived due to the property of the operator 
of P t . Since the operator P e is defined on L 2 (0, T; H k (D e ), k = 0, 1, for / Iff. 6}) 
we define 

P e ii e = A e P £ u £ + f e + g e W, onDx (0, T). 
By the property of P e and the estimates of u e , it is easy to see that 

P t u e = {PeU e ), inD e x (0,T) 

and 

E|P e 'U e |l / 2(o i T;ii"-l) < E|-u e | L 2( 0)T;He -i). 

4. Effective macroscopic model 

We now derive the effective macroscopic model for the original model (12.11) . 
Let u e £ L 2 (0, T; H e ) be the solution of system (jO!- (j2^jl . Then by the esti- 
mates in Theorem 13.11 Remark 13.51 and the Chebyshev inequality [121 H3] , it is 
clear that for any 5 > there is a bounded set K$ C G with spaces X, y and 
Z in Lemma 12. II (and in the paragraph immediately before it) are replaced by 
Hq(D), H and H~ l (D) respectively, such that 

¥{P e u e £ K s } >l-5. 

Thus Ks is compact in L 2 (0, T; H) by Lemma I2TT1 Then {C(P t u t )} t is tight in 
L 2 (0,T; H). The Prokhorov Theorem and the Skorohod embedding theorem 
([T2"]) assure that for any sequence {e^} with €j — >• as j — > cxo, there exists 
a subsequence {ej(fc)}, random variables {u e . (k) } C L 2 (0,T; H e . (k) ) and w £ 

L 2 (0,T; if) defined on a new probability space (O, J 7 , P), such that 
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and 

P e . (k) u em ^u in L 2 (0,T;H) as k -> oo, 

for almost all u G Vt. Moreover P ej . k ,u e ., k) solves system (I2.ip - (12.4p with W 

replaced by Wiener process Wk defined on probability space (fi,.F, P) with 
same distribution as W. The limit u is unique; see jl], p. 333. In the following, 
we will determine the limiting equation (homogenized effective equation) that 
u satisfies and the limiting equation is independent of e. After this is done we 
see that C(u e ) weakly converges to C(u) as e | 0. 

We always assume the following 

fe^f, weakly in L 2 (0,T;H) (4.1) 

and 

gl^g\ weakly in L 2 (0,T;H). (4.2) 
Define a new probability space Fs, P$) as 

VL$ = {uj E fi : u e {uS) E K$}, 

F s = {Fnn s :FeF} 

and 

w-W- 

Denote by E5 the expectation operator with respect to Fg. 

Now we restrict the system on the probability space (Q5, P5). In the 
following discussion we aim at obtaining L 2 (Q$) convergence for any 5 > 
which means the convergence in probability [31 [H] . 

From the estimates (13.51) . (I3.6p . Remark I3~51 and the compact injection G C 
L 2 (0,T; H), there exists a subsequence of u e in Kg, still denoted by w e , such 
that for a fixed u E Q5 

P t u e u weakly* in L°°(0, T; H) (4.3) 
P e u e u weakly in L 2 (0, T; iJ 1 ) (4.4) 
strongly in L 2 (0,T;H) (4.5) 
P t u e ^ii weakly in L 2 (0,T; H' 1 ). (4.6) 



Define 
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which satisfies 

- div^ e = f e + g e W - u t in D e x (0, T) (4.7) 

£ e . n = o ondS e x (0,T). (4.8) 

By the hypothesis of and the fact that {u e )e being bounded in L 2 (0, T; Hq), 
we have 

£ £ ->> £ weafc/y m L 2 (0, T; © n #). (4.9) 

We make use of Tartar's method of oscillating test functions to determine the 
limiting equation [9]. 
Note that 

(*T p pT p ^ pT p 

£ e -Vvifdxdt = f e vtpdxdt + y^ / / g\vdx(pdWi{t) + 

Jo J d i=1 Jo Jd 

T p 

P e UeXD^vdxdt (4-10) 



JD 



o Jd 

for all v G Hl(D) and <p E V(0, T). We pass to the limit in (OOj) as e -> 0. 
Due to the fact 

P e u e ^u strongly in L 2 {0,T;H), (4.11) 
Xd c ^$ weakly* in L°°(D) (4.12) 

and the estimate 

OO n-J 1 „ „ OO 



~i|2 I |2 

e I L 2 (0,T;H) I ^ I L 2 (0,T;H) > 



1 i=l ^ i=l 

by the assumption (14.21) . we see that 

V / / glvdxipdWi{t) -»• V / / j'wd^d^fi), m L 2 (ft). 

. -, Jo Jd ,_i Jo Jd 



t 'o Jd i=1 jo 

Thus letting e — > in (I4.10p and since L 2 (fls) is a subspace of L 2 (f2) one 
finds that in L 2 (Qs) 

pT p pT p °° pT p 

' / £ • VvLpdxdt = fvLpdxdt + / / g l vdx(pdWi(t) 

o Jd Jo Jd i=1 Jo Jd 

+ [ [ ®u(pvdxdt. (4.13) 
Jo Jd 

Hence 

- div £(x, t) = f(x, t) + g(x, t)W - $u in D x (0, T). (4.14) 
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In the following we identify the limit £. We follow the approach of deter- 
ministic case for the elliptic problem with homogeneous Neumann boundary 
condition [9]. 

For any A G M n , let w\ be the solution of 

n a n a 

-x;|-(|><fH^* Ms) 

w\ — \ ■ y Y — periodic (4.16) 
<9^a 



on 95 (4.17) 



and define 



fx 

w{ = e(Qw A )(- 



where Q is in Lemma I3.2L Then we have [9] , 

w\ A • x weakly in H l {D), (4.18) 
Vw\ A wea£% in ® n L 2 (D). (4.19) 

Now we define 

v i=i ayi J 

and (rjt)(x) = {r]^{x/e))j = A\Vw\. Then 

_ di v ^ A = in D (4.20) 

and due to (14181) and 11051) 

r? A My(r] A ) wea£% m L 2 (D). (4.21) 

It is easy to see that J\4y(t) X ) = B t X with 5* = {(3ji) a constant matrix which 
is determined in the appendix. 

Using test function (pvw\ with (p G X>(0, T), u G T>(D) in (14.101) and multi- 
plying both sides of (14.201) with ipvP e u e , one has 



/ / £e • Vvipw e x dxdt + / £ e • Vw\v<pdxdt 
Jo Jd Jo Jd, 



rf e ■ Vv(fP e u e dxdt - \ rf e ■ V (P e u e )v(pdxdt 
o Jd Jo Jd 



00 /»T /» rT I* 

f e ipvw\dxdt+y / / glvw t x dx(pdW i (t)+ / / P fL u fL XD,, ( P vw \dxdt. 
i=1 io Jo Jd 



T p 00 /-T 

o 
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Then by the definition of £ £ , and the assumptions ( 14. ip . ( 14. 2p . using the 
convergence (gU), (|4TTT]h l|OZ]l . (l4TTgl) . (Q5]l and (ICTj) . we have in L 2 (n 5 ) 

/ / £ ' Vf<£>A " xdxdt — I I B t X- Vvipudxdt 
Jo Jd Jo Jd 

pT p pT p pT p 

= / fpv ^ " xdxdt + / / (7 i t'A-X(ia;(^(iWj(t)+ / / fluvLpX- xdxdt. 
Jo Jd i=1 Jo Jd 



o 



That is 



/ / £ ' V(f A ■ x)cpdxdt — I / £ • XvLpdxdt — I I B l X ■ Vv^pudxdt 

Jo Jd Jo Jd Jo Jd 

/ / fipv X ■ xdxdt + 2^ / / g l vX-xdx(pdWi{t) + 
Jo Jd i=l Jo Jd 



oo r T r r T f 

I -QuvtpX- xdxdt. 
o Jd 



Then by using (14.131) with the test function replaced by v X ■ x^p one has 



T r rT 

t • 



£ • Xvipdxdt = I B X ■ Vmpvdxdt 
o Jd Jo Jd 

which yields 

£ • A = B l X ■ Vu = BVu ■ A. 

Then 

f = BVu 

since A is arbitrary. Then u satisfies the following equation 

Mu = (div(BVu) + f)dt + gdW(t). (4.22) 

Suppose 

u° e u°, weakly in H. (4.23) 

We now determine the initial value by suitable test-functions. In fact, taking 
v e V(D) and <p e V([0,T]) with <p(T) = we have 

T p pT p 00 pT p 

£ e ■ Vvipdxdt — I I feVtpdxdt + / / glvdx<pdWi(t) — 
Jo Jd -n Jo Jd 



o Jd 



8=0 

,~.o, 



u € vipdxdt + / u e ip{Q)vdx. 
o Jd Jd 
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Now let e — > 0, noticing that 

/ u t vipdxdt = XDcPeUeVtpdxdt — > / / -duwpdxdt 

a J d Jo Jd Jo Jd 

/ fiiivLpdxdt + / $u(0)(p(0)vdx 
o Jd Jd 



by (14.141) . we have 



u° 



Here one should notice that the above result is in the sense of L 2 (Qs)- Then 
the above analysis yields the following results 

limE 5 |P £ u e - u\ 2 L2{0T . H) = (4.24) 

and 

limE 5 / (A e P e u e - BVu)v(pdxdt = (4.25) 
Jo Jd 

for any u G and G X>([0,T]). 

Now we are in the position to give the homogenized effective equation in the 
following theorem. 

Theorem 4.1. (Effective macroscopic model) For any T > 0, assume 
that ffljfy and j4.23\j hold. Let u t be the solution of If2l\) -$T4\). Then the 



distribution C(P e u e ) converges weakly to fi in the space of probability measures 
on L 2 (0, T; H) as e { ; with fi being the distribution ofu, which is the solution 
of the following homogenized effective equation 

Mu = (div(BVu) + f)dt + gdW{t) in D x (0, T), (4.26) 

u = ondD x (0,T), (4.27) 
u° 

u(x,Q) = — inD, (4.28) 

17 

where B = (fy) is determined by fl7.^| ) in Appendix at the end of this paper. 

Remark 4.2. This theorem implies that the macroscopic model \4-26$ is an 
effective approximation for the microscopic model \2. on any finite time 
interval < t < T , in the sense of probability distribution. In other words, 
if we intend to numerically simulate the microscopic model up to finite time, 
we could use the macroscopic model as an approximation when e is sufficiently 
small. 
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Remark 4.3. Due to the appearance of the stochastic integral term (see M.lOj) ), 



this theorem on weak convergence of probability measures does not follow di- 
rectly from the deterministic homogenization results and the mild formulation 

MM- 



Remark 4.4. The stochastic PDE 1^4- 261) is defined on the homogenized do- 
main D. By the analysis in [12], for any fixed T > 0, the macroscopic syste 
n^-njm is well-posed, as long as f G L 2 {0,T;H) and g G L 2 (0,T;£^). 



em 



Proof. Noticing the arbitrariness of 5, this is direct result of the analysis of the 
first part in this section by the Skorohod theorem and the L 2 (Qs) convergence 
of P e u e on (fi tf ,^,Pj). □ 

We finish this section by the following remark. 

Remark 4.5. Note that there are several papers on effective dynamics for 
partial differential equations with random coefficients (so called random PDEs; 
not stochastic PDEs); see [19j [26j [32] and reference therein. In [HI [26], a 
random partial differential equation is obtained as the homogenized effective 
equation for a random system with fast or small scales on both time or spatial 
variable. And the distribution of solution of heterogeneous system converges 
weakly to that of homogenized equation. However in [32], the effective equation 
is obtained as an averaged deterministic equation for a random system with 
small scale just on time. And the fluctuation of the solution of the random 
equation around the solution of the averaged equation converges to a generalized 
Ornstein- Uhlenbeck process in distribution. In the present paper, the original 
microscopic model is a stochastic PDE (i.e., PDE with white noise) and the 
effective macroscopic equation is still a stochastic partial differential equation. 

5. Long time effectivity of the macroscopic model 

In this section we consider the long time effectivity of the homogenized sys- 
tem (I4.26P in the autonomous case. It is proved in section H] that for fixed 
T > the macroscopic behavior of the microscopic system (12. ip ~ (12.41) can be 
approximated by the macroscopic model (I4.26P in the sense of probability dis- 
tribution. In fact we can show the long time approximation. More specifically, 
we now prove that in the sense of distribution, all solutions of (I2.1l) - (l2.4p con- 
verge to the unique stationary solution of (I4.26P as T — > oo and e — > 0, under 
the assumption that f e G H e and g\ G V e are independent of time t and 

oo 

J2\VA e 9l(x)\l nHe <C*. (5.1) 

i=l 

Here C* is a positive constant independent of e. 
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By the above assumptions, the property of ctij and fy, a standard argu- 
ment (see [IB] , Section 6) yields that the system (13. ip and (I4.26P have unique 
stationary solutions u*(x, t) and u*(x, t) defined for t > 0. We denote by /i* 
and /i* the distributions of P e u* and w* in the space H, respectively. Then if 
E|w°| 2 < oo and E|w°| 2 < oo, 



hdfi e (t) 



H 



H 



<C(u° t )e-^, t>0, 



(5.2) 



hdfj.it)- I hdfi* <C(M u )e" 7 *, t>0 



(5.3) 



for some constant 7 > and any h : H — > R 1 with sup \h\ < 1 and Lip(/i) < 1. 
Here fi e (t) = C(P e u e (t,u° e )), n{t) = and C(u° e ) and C(u°) are 

positive constants depending only on the initial value vP e and u° respectively. 
The above convergence also yields that /i e (t) and nit) weakly converges to fi* 
and fi* respectively, as t — > 00. 

We will give some additional a priori estimates which is uniform with respect 
to e to ensure the tightness of the stationary distributions. For Banach space 
U and p > 1, we define W 1,p (0, T; U) as the space of functions h e L p (0, T; £7) 
such that 



W 1 -P(0,T;U) 



\h\ p 4- 



eft 



LP(0,T;f/) 



< OO. 



And for any a e (0,1), define W a ' p (0, T; ?7) as the space of function /i e 
L p (0,T;£f) such that 



W>P(0,T;C7") 



f - M+ap 



For p G (0, 1), we denote by C p (0,T; U) the space of functions h : [0,T] 
that are Holder continuous with exponent p. 



X 



In the following part of this section we always assume that f t and g\ are 
independent of time t with (15. ip hold. And for T > denote by u* T (respec- 
tively, Uy) the distribution of stationary process P £ u*(-) (respectively, w*(-)) in 
the space L 2 (0,T; H 1 ). Then we have the following result. 

Lemma 5.1. For any T > £/ie family u* T is tight in the space L 2 (0, T; H 2 ^ i ) 
with 1 > 0. 

Proof. Since w* is stationary, by Q3.8j) . we see that 



E|m*|l2( 0jT;jH -2) < C T . 



(5.4) 
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Now represent u* in the form 

u*{t)=u*{0) + f A e u*(s)ds+ [ f e {x)ds+ [ g e {x)dW{s). 
Jo Jo Jo 

Also by the stationarity of u* and (13.81) we obtain 
E 



A e P e u*(s)ds+ / f e (x)ds <C T . (5.5) 

Jo W^(0,T;H) 

Let M t (s,t) = p g e {x)dW \s) . By Lemma 7.2 of [12] and Holder inequality, we 
derive that 

E|M e (M)& ^c^j* \V A Mx)\ 2 c Qdr) 2 < K(t - s) j* |V aM^^t 

< KC* 2 \t-s\ 2 

for t e [s, T], where K is a positive constant independent of e, s and t. Then 

E f \M e (0,t)\ 4 v dt<C T (5.6) 
Jo 

and 

Combining fl5.4l) - (l5.7p . and the compact embedding of 

L 2 (0, T; F 2 ) n W /1 ' 2 (0, T; H) C L 2 (0, T; if 2 "') 

and 

L 2 (0, T; if 2 ) n W aA {0, T- H l ) C L 2 (0, T; /J 2_t ) 
we obtain the tightness of u* T . This completes the proof. □ 

The above lemma directly yields the following result 

Corollary 5.2. The family {/i*} is tight in the space H 1 . 

By Lemma 15.11 for any fixed T > 0, the Skorohod embedding theorem 
asserts that for any sequence {e n } n with e n — > as n — > oo, there is sub- 
sequence {e„(,t)}fc, a new probability space (fl, J 7 , P) and random variables 
^ G L 2 (0, T; K), m* e L 2 (0, T; H 1 ) such that 

and 

«^, fc) — ► «*, in L 2 (0,T;iJ 1 ) as /c — > oo. 
Moreover w* n{ft) (respectively, «*) is the unique stationary solution of equation 
(13.11) (respectively, (I4.26P ) with W replaced by Wk (respectively, W). Here Wk 
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and W are some Wiener processes defined on (O, J 7 , P) with same distribution 
as W. Then by the analysis of Section H] and the uniqueness of the invariant 
measure, we have 

u* T — u^, as e — > 

for any T > 0. 

To show the long time effectivity, let u e (t), t > 0, be a weak solution of 
system ( I2.ip - (I2.4I) and define «*(•) = u e (t + -) which is in the space Lf oc (R + ; V e ) 
by Theorem 13.11 Then by (15. 2p 

C(P e ul(-)) ^ C(P e u* e (-)), t^oo 

in the space of probability measures on L^ oc (R + ; H 1 ). Having the above anal- 
ysis we draw the following result which implies the long time effectivity of the 
homogenized effective equation (14.261) . 

Theorem 5.3. (Long time effectivity of macroscopic model) 

Assume that f e G H t and g\ G V e are independent of time t with 115.1}) hold, 
and further assume that fl^.l| ) and Iji4-S\ ) hold in H . Denote by u e (t), t > 0, the 
solution of I12.1\) - (K4\ ) and u* the unique stationary solution of \4.20\j . Then 

Iimlim£(P e tt*(.)) =£(«*(•)), (5-8) 

ej.0 t^oo 

where the limits are understood in the sense of weak convergence of Borel prob- 
ability measures in the space Lf oc (M. + ; H 1 ). That is, the solution of 112. 

converges to the stationary solution of ^4-26\) in probability distribution 
as t — > oo and e — > 0. 



Remark 5.4. This theorem implies that the macroscopic model H4-26] ) is an ef- 
fective approximation for the microscopic model (OOP; on very long time scale. 
In other words, if we intend to numerically simulate the long time behavior 
of the microscopic model, we could just simulate the macroscopic model as an 
approximation when e is sufficiently small. 



6. Effectivity in energy convergence 

In the last two sections, we have considered finite time and long time effectiv- 
ity of the macroscopic model (I4.26p . in the sense of convergence in probability 
distribution. In this section we focus on the finite time effectivity of the macro- 
scopic model (I4.26p . but in the sense of convergence in energy. Namely, we 
show that the solution of the microscopic model (12.11) or (13. ip . converges to 
the solution of the macroscopic model (I4.26p . in an energy norm. 



and 
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Let u e be a weak solution of ( 13. ip and u be a weak solution of ( I4.26p . We 
introduce the following energy functionals: 

£ e «>(*) = \^\u e \ 2 H + E / / X De^eV(P e w e (a; ) r))V(P e « e (x ) r))dxdT (6.1) 

£°(u)(t) = ^E\u\ 2 H + E [ [ BVu(x,T)Vu{x,T)dxdT. (6.2) 
2 Jo Jd 

By the Ito formula, it is clear that 

£ £ {u e )(t) = ]-E\u° e \ 2 H + E [ [ f e {x,T)u e {x,T)dxdT + ]-E [ |# e (x, r)| 2 Q dr 
^ Jo Jd ^ Jo 2 

and 

= J E k% + E / [ f{x,r)u{x,T)dxdr + lE [ \g(x,r)\ 2 Q dr. 
2 Jo Jd * Jo 2 

Then we have the following result on effectivity of the macroscopic model 

in the sense of convergence in energy. 

Theorem 6.1. (Effectivity in energy convergence) 
Assume that UTlty and fiJJty hold. If 

— > u°, strongly in H, as e — > 0, 

then 

£ e (u t )^£°(u) in C([0,T)),ase^0. 

Proof. By the analysis of section HJ for any 5 > 0, u t — * u strongly in 
L 2 (0, T; H) on f^, then by the arbitrariness of 5, it is easy to see that 

E / / fe(x,T)u e (x,T)dxdr — ► E / / f(x,T)u(x,r)dxdT, fort£[0,T]. 
Jo Jd Jo Jd 

Then by g e — ^ (7 weakly in L 2 (0, i; £2 )> we have 

£ e (u e )(t) ->£°(«)(t) for any t E [0,T]. (6.3) 

We now only need to show that {£ e {u e ){t)} t is equicontinuous, as then the 
Ascoli-Arzela's theorem [H] will imply the result in the theorem. 
In fact, given any t G [0,T], and h > small enough, we have 

\£ e (u e )(t + h)-£ e (u e )(t)\ 

/t+h r _ rt+h 
J f e (x } T)u e (x,T)dxdr \9e{ x , T )? c QdT 

< Ejl/el^^Tjfl) y |M £ (x,r)|^rfa;rfr| + Ey |^ e (x, r)|Jgdr. 
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Noting that u e G L 2 (0, T; H) a.s. and (12.71) . we have 

\£ £ (u t )(t + h) - £ e {u e ){t)\ ->0, ash^O, 

uniformly on e, which means the equi-continuity of the family {S e (u e )} e . This 
completes the proof. □ 

7. Appendix: The homogenized matrix 

In this Appendix, we give the explicit expression of the homogenized matrix 
B; for more details see [9]. Let x\ i — lj • • • > n be the solutions of 

-Ef U^zM) = o y (7.i) 



V °fc< a w ' =0on9S (7.2) 

X % is K — periodic. (7.3) 

It is easy to calculate that x l — ~ w e i + e« with {e;}™ =1 the canonical basis of 
R n . Then 

& = V\Jy^ ak °^ dV = V\Jv ^ ~ M X £ ^ 

Moreover the operator i? = satisfies the uniform ellipticity condition: 

there is a constant 6 > such that 

2 > 6 E & M £ = (6, • • • , £n) e K". 

ij'=l i=l 
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